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Abstract. We consider a linear transport equation on the edges of a network with time- 
varying coefficients. Using methods for non-autonomous abstract Cauchy problems, we ob- 
tain well-posedness of the problem and describe the asymptotic profile of the solutions under 
certain natural conditions on the network. We further apply our theory to a model used for 
air traffic flow management. 



1. Introduction 

Dynamical processes taking place in networks have been of enormous interest in recent 
years and have various applications for real life phenomena. We are interested in transport 
processes or flows in networks. Methods from the theory of operator semigroups to treat such 
processes were first used in |KS05| for a finite network where a simple transport equation 

d . , e? 
—u(x,t) = Cj —u(x,t) 

was considered on the edges together with boundary conditions of Kirchhoff-type in the 
vertices. These methods were further applied to various generalizations of this problem in 
finite |Sik051 IMS071 IRad08j or even infinite networks |Dor08L IDKS09| . The authors obtain 
well-posedness and describe the asymptotic behavior of the solutions. Further, [EKNS08, 
EKKNS10J studied control probl ems for flows in networks. See also [DKNR10J for a survey 
of the semigroup approach to transport processes in networks. 

The processes in all mentioned works are autonomous, i.e. the differential operators gov- 
erning the processes do not change in time. Motivated^ by applications to air traffic flow 
management (see Section [5]), we now study non-autonomous processes. Precisely, we are in- 
terested in transport processes where the boundary conditions in the vertices vary in time. 
This yields differential operators with varying domains and the corresponding Cauchy prob- 
lems become non-autonomous. Solutions to such problems are described by evolution families 
instead of semigroups, see |Bayl2a[ IEN00| |Nag95[ INN02] . 



In the following we first define the time-depending network with a transport process in it. 
Our main tool to study such non-autonomous processes is the theory of difference evolution 
equations as developed in |Bayl2a |Bayl2b| which we briefly describe in [3j The main results 
are contained in Sections [4] and [5 where we treat two different flow processes in a network 
and prove well-posedness of both of these problems. Assuming periodic boundary conditions, 
we obtain asymptotically periodic behavior of the solutions. The period is given in terms of 
the (time-depending) network structure. 

The last author would like to acknowledge the support of the German Academic Exchange Service (DAAD) 
during her stay at the University of Tubingen. 

"*"We are grateful to Benedetto Piccoli for drawing our attention to this problems. 
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2. Preliminaries 

2.1. Time-depending networks. The network is modeled by a finite directed graph G 
consisting of n vertices Vi, . . . ,v n and m directed edges (arcs) ei, . . . ,e m . We equip every 
edge ej with time- varying weight uiij{t) > such that 

m 

(1) y^ufy-ft) = 1 for every t G M + and every i 

i=i 

(here i numbers either vertices or edges — it depends on the concrete problem and we will 
specify it later on). The graph structure is described by the outgoing incidence matrix & = 

with 



1, ifvj— 
0, otherwise, 



and the incoming incidence matrix & + = (<Ptj) with 



nxm 
e, 



3+ 



if 

otherwise. 



Instead of using incidence matrices, it is sometimes more convenient to use adjacency matrices. 
Here, we use the (transposed) adjacency matrix of the line graph B = (bij) mxm with entries 



if — > v — 
otherwise. 



A directed graph is called strongly connected if for any pair of distinct vertices Vj , Vj there 
is a directed path in the graph going from Vj to vj and vice versa. This property can be 
characterized by irreducibility of the usual vertex adjacency matrix (see e.g. |Min88l Theorem 
IV.3.2]), but also by our adjacency matrix of the line graph. 

Lemma 1. |Dor08| Proposition 4.9] A directed graph is strongly connected if and only if the 
matrix IB is irreducible. 

2.2. Transport processes. In order to model a transport process on the edges, we normalize 
the edges as ej = [0, 1] and parameterize them contrary to the direction of the flow, i.e., the 
material flows from 1 to 0. We consider some finite mass distributed on the edges of the 
network and denote by Uj(x, t) its density at position x G [0, 1] of the edge ej and at time t, 
hence Uj : [0, 1] x R — > R, j = 1, . . . m. 

Our basic assumptions on the process are the following. 

(1) On each edge ej we describe the transport process by 

— Uj (x,t) = — . 
at ox 

(2) The initial distribution of the mass on the edges ej at time s G R is given by 



-Uj (x, t) = -jr-Uj (x, t) , x G (0, 1), t > s. 



Uj{x,s) = fj(x), x G (0, 1). 
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(3) No mass is gained or lost during the process. In particular, no absorption takes place 
along the edges, and in each node Vj we have a Kirchhoff law 

m m 
3=1 k=l 

(4) In each vertex Vj the incoming material is distributed into the outgoing edges 
according to the time- varying weights uJij(t) > so that holds. 

By choosing two different ways to assign the weights tOij(t) to the edges, we will in Sections 
[4] and [5] obtain two different flow processes in the network. In the first case we will assume 
that the material is collected in the vertex and is then redistributed according to the weights. 
In the second case we want to keep track of the origin of the material and hence the weights 
will give the proportions of the material that flows from one edge into another one. 



3. NON-AUTONOMOUS DIFFERENCE EQUATIONS 

To tackle our transport problem in time-depending networks we will use the theory of 
positive evolution families corresponding to a class of non-autonomous difference equations 
developed in |Bayl2a| and |Bayl2b| . We explain the terminology and state the results needed 
below. 

Choose the Banach space X = L 1 ([0, l],C m ) as the state space of the system. For a family 
of matrices (-B(i)) tgK C M m (C) we define difference operators A(t) : D (A(t)) — > X by 

(2) D (A(t)) := {/ G W 1 ' 1 ([0, 1], C m ) | /(l) = B(*)/(0)} and A(t)f := f 

for / G D (A(t)) and t G K. The non-autonomous abstract Cauchy problem corresponding to 
the operators (A(t), D (A(t))) is of the form 



(nACP) 



u(t) = A(t)u(t), t>s, 
[ u(s) = fs e X. 

A classical solution to the (nACP) is a differentiable function u G C 1 ([s,oo),X) such that 
u(t) £ D (A(t)) for every t > s and u satisfies (nACP). Furthermore, we say that (nACP) is 
well-posed if there exists a unique evolution family (U(t, s)) t>s such that that the regularity 
subspaces 

Y s := {/ G X | [s,oo) 3 1 1 ^ U(t,s)f is a classical solution to (nACP)} 

are dense in X for every s G K. For the definition of the evolution family see |Bayl2a| 
or [ ENOOl Definition VI.9.2]. We also recommend |JNag95| , |JNM02j . or |Paz83j for further 
information on evolution families and their relation to non-autonomous Cauchy problems. 

Since the domains D (A(t)) are time-dependent and do not contain a common core, none 
of the usual well-posedness results is applicable in our case. We will use the following results 
from |Bayl2a| instead. 

Proposition 2. |Bayl2a Theorem 2] Let the mapping t \— > B(t) be uniformly bounded and 
absolutely continuous. Then the (nACP) associated to the operators (A(t) , D (A(t))) given by 
^ is well-posed. 

In |Bayl2a| even an explicit formula for the corresponding evolution family is given. We 
state here this formula in a special case. 
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Proposition 3. |Bayl2a[ (8)] Let the mapping t \— > B(t) be uniformly bounded, absolutely 
continuous and 1-periodic, i.e. B{t + 1) = B(t) for every i £ 1. Then the unique classical 
solution to {nACP) is given by u(t) = U(t, s)f s , where 

(3) U(t,s)f = B k {t + x)f{x + t- s-k), 

for f £ X, x G [0,1], k < x + t - s < k + 1 and k G N . 

In order to study the asymptotic behavior of the solutions, some more regularity assump- 
tions are needed. Denote the unit circle by T := {z £ C | \z\ = 1}. 

Proposition 4. Let t \— > B(t) be a continuous, 1-periodic mapping and let B(t) be a stochastic 
and irreducible matrix for every t£l. Then there is a family of projections {P(s) \ s G M} 
in C(X), commuting with the evolution family (U(t, s)) t>s and decomposing the space X as 

X = X n (s) X s {s) := P(s)X kerP(s) 

such that the following properties hold. 

(i) The subspaces X-ji(s) and Xs(s) are (U(t, s)) t>s -invariant for every sgM. 

(ii) (Us(t, s)) t>s := (U(t, s)\x s ( s )) t> is uniformly exponentially stable, i.e. there exist 
C > 1 and uj > such that 

\\U S (t,s)\\ <Ce-^- s \ t>s. 

(Hi) (Ufc(t, s)) t>s := (U(t,s)\x n ( s )) t>s can be extended to an invertible evolution family 
(U-ji(t, s))u s \ e ^2 which is positive and periodic in evolution, i.e. Uh(s + t, s) = Ix n (s) 
for every s G M, with the period 

t = icm{|5(t) nr| | t G [0,1]}. 
(iv) For t as above there exists a r-periodic positive group (T(t)) 4gR such that 

\\U(t,s)-T(t-s)P(s)\\ t ^0 

for every s£R. 

Proof. The m x m matrices B(t) are all stochastic and irreducible, therefore by Perron- 
Frobenius theory (see |Sch74| Theorem 1.6.5]) the peripheral spectrum a (B(t)) n T, for every 
t G M, is a finite group consisting of (at most m) roots of unity which are all first order poles 
of the resolvent. Hence the union 

|J {a(B(t))nr} 

te[o,i] 

is a finite discrete set and the least common multiple 1cm in (iii) is well defined. The stochas- 
ticity of the matrices B(t) also implies that the evolution family (U(t, s)) t>s given in Q 
consists of positive contractions. Combining Proposition 6, Definitions 7 and 8, and Theorem 
9 from |Bayl2a| we now obtain the decomposition of the space X with properties (i) and (ii). 
Finally, (iii) follows from )B ayl2a[ Theorem 25] and (iv) from |Bayl2a Theorem 26]. □ 
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4. Flows in nonautonomous networks 



Consider now a finite weighted network G as in Section 2.1 with incidence matrices $~ 
and <1> + . The time-dependent weights 0Jij(t) > in every vertex Vj give the proportions of the 
incoming material to be distributed into the outgoing edges at time t, where 

Wy(t) = 0if 0y =0. 

This condition reflects the fact that the edges of our network are fixed and the flow takes 
place only on the edges of the network. Note however that it might happen that no material 
is sent from the vertex Vj into the edge ej at time to for some to, meaning that Wij(io) = even 
if ^ 0. We store the weights in the time-dependent weighted outgoing incidence matrix 
= (<?W(*))nxm defined as 

0- ..(f) : = J w y(*)' ifv *^' 
w ' tJ lo, otherwise. 

The m x m time-dependent weighted adjacency matrix of the line graph M w (t) is obtained by 
(4) M w (t):={<S>-(t)) T $ + - 

Note that the nonzero entries of M w (t) are in one-to-one correspondence with the nonzero 
entries of the unweighted adjacency matrix IB of the line graph. We assume that there is no 
absorption in the vertices, hence ([I]) holds for every i G {1, . . . , n} and all t G M+, and the 
matrices M w (t) are all column-stochastic. 

By Gt we will denote the network at time t obtained from the adjacency matrix M w (t). 
This means that Gt C G where the edges of G with no inflow at time t are deleted. 

Under these assumptions we study the following transport process in G. 

§tUj(x,t) = ?fcUj(x,t), x G (0,1), t > s, 
Uj (s,Q) = £00,5 €(0,1), (IC) 

br jUj (l,t) = ^(tJEH^C.*):^ ( nBC ) 

for i = 1, . . . , n and j = 1, . . . , m. It is of the form given in Section |2.2| Note that the non- 
autonomous boundary conditions (nBC) together with ([!]) imply the Kirchhoff law 2.2 (3j). 



(nF) < 



In order to use the results from Section[3j we now take the Banach space X = L 1 ([0, 1], C r 
and the difference operators A^(t) on X associated to the family of matrices (B u ,(t)) tgK as 
defined in ([2]), hence 

An(t)f := /' with domain D (A M (t)) := {/ G W 1,1 ([0, 1], C m ) \ /(l) = B w (t)f(0)} . 

Proposition 5. The non- autonomous abstract Cauchy problem (nACP) corresponding to 
(A^,D (Aa(t))) is equivalent to the transport process in the time- depending network (nF). 

Proof. We only need to observe that the non-autonomous boundary conditions (nBC) of the 
problem (nF) are hidden in the domain D (^4b(^)) : 

g G W 1 ' 1 ([0, 1], C m ) satisfies (nBC) <J=> g(l) = B w (t)g(0), 

similarly as in |Dor08} Prop. 3.1]. □ 

The well-posedness of the analogous problem in the autonomous case (even for infinite 
networks) accompanied with an explicit formula for the solution was shown in [Dor08, Prop. 
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3.3]. By Propositions [2] and [3] we can prove a well-posedness result of our non-autonomous 
flow problem (nF) for a special class of periodic time-dependent networks. 

Corollary 6. If the mappings t >— > ooij(t) are absolutely continuous, then (nF) is well-posed. 
In particular, if in addition, the mappings t 1— > M w (t) are 1-periodic, i.e. M w (t + 1) = M w (t) 
for all t G M, then the unique classical solution to the non- autonomous flow problem (nF) in 
G is obtained by the evolution family in |3]) as 

(5) u(t) = U(t, s)f s =M k w (t + x)f s (x + t-s-k), 

where f s is given by the initial conditions (IC), t > s, x G [0, 1], k < x + t — s < k + 1, and 
/c € No. 

Note that the period is assumed to be 1 only as a matter of convenience. We could take 
any natural number and modify the above formula appropriately. 

Using Proposition [4] we obtain an asymptotically periodic behavior of our non-autonomous 
flows. 

Theorem 7. Let the network Gt be strongly connected for every t G 1 and let the mapping 1 1— > 
B UI (t) be continuous and 1-periodic. Then the flow evolution family ^ converges uniformly 
to a periodic positive group in the sense of Proposition^ Its period r can be computed as 

t = lcm {gcd{/ I ej 1 , . . . &j l form a cycle in Gt} \ t € R} . 

Proof. By ([I]), the matrices M w (t) are all column-stochastic. Since the graphs Gt are all 
strongly connected, the matrices B„,(i) are all irreducible by Lemma [TJ Hence we can apply 
Proposition [4| For the expression for the period r note that \a (M> w (t)) n T| equals the index 
of imprimitivity of the matrix M w (t) (see |Min88, Definition III. 1.1]), which is the same as 
the greatest common divisor of all cycle lengths in the network Gt, cf. |Min881 Theorem 
IV.3.3]. □ 

Remark 8. If we assume that the nonzero weights remain strictly positive in time, i.e. 

0- ^ =>Uij(t) ^ forallt€R, 

then the index of imprimitivity ofM w (t) does not change in time either, and the period r can 
be computed simply as 

t = gcd{Z I ej 1 , . . . &j l form a cycle in G}. 

This means that in the case of non- disappearing edges in the network, the strictly positive 
weights do not have any impact on the period and the asymptotic behavior remains the same 
as in the autonomous case (see [KS05, Corollary 4.7],). 

4.1. Examples. 

(1) We consider the family of networks Gt depicted in Figure [l] The weights on the 
edges ei, e2, e3, e6 are constant, the weights on the edges e± and es vary in time. 
Each edge has a nonzero flow of material on it at every time. The corresponding 
adjacency matrix M w (t) is 
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Figure 1. The flow on this network is asymptotically periodic with period 1. 
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The mapping 1 1— > M w (t) is absolutely continuous, so by Corollary [6j the correspond- 
ing non-autonomous network flow problem (nF) is well-posed. The matrices ~B w (t) are 
all column stochastic, 1-periodic, and, since the graphs Gt are strongly connected for 
all t E R, irreducible. We can therefore apply Theorem[7j Note that we are in the case 
of Remark [8] since each edge carries a nonzero flow of material at every time t 6 M. 
For every t, the graph Gt contains a cycle of length 3 and a cycle of length 4. Hence, 
the flow evolution family converges uniformly to a periodic positive group with period 
r = gcd{3,4} = 1 in the sense of Proposition [4j meaning that the flow on this graph 
asymptotically behaves periodically with period r = 1. 
(2) The networks in Figure[2]all belong to the family of networks Gt given by the adjaceny 
matrix 
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Note that the weights on edges eg and eio do not vary in time, but those on the 
remaining edges do. This is indicated in Figure [2] (i) . Hence, some edges do not carry 
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Figure 2. The three possible states of the family of networks Gt in Example 2: 
Either, all edges of the network carry flow (Figure (i)), or there is no flow on 
the edges belonging to the outer cycle of length 4 (Figure (ii)) or to the inner 
cycle of length 4 (Figure (hi)). 



any flow for some t 6 R: For t £ Z, the weights on the edges es, e&, ej, es are zero 
(Figure I2J (ii)), whereas for t G ^Z \ {0}, there is no flow on the edges ei, &2, 63, e4 
(Figure pj(iii)). Since eg, eio form a cylce of length 2, the greatest common divisor of 
all cycle lengths in the network is equal to 2 for all times f G K, and hence also the 
least common multiple appearing in Theorem [7] is equal to 2. 

Again, we obtain well-posedness of the corresponding flow problem (nF) by Corol- 
lary [6] since the mapping t 1— > M w (t) is continuously differentiable. Note that the 
graphs Gt are strongly connected for all t £ M and fulfill all assumptions required by 
Theorem [7j implying that the flow evolution family converges uniformly to a periodic 
positive group with period r = 2. 
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5. Application to air traffic flow management 



The following application of our results is motivated by the (real world) regulation of air 
traffic, called Air Traffic Flow Management. Its goal is to optimize air traffic flow, i.e., limiting 
the density of aircraft in certain regions of airspace as well as operating efficient routes subject 
to weather constraints. These tasks are currently prescribed by playbooks established over 
time and based on controller experience. However, one of the aims consists in providing a 
mathematical model of air traffic flow allowing to apply mathematical control techniques. 

5.1. Modelling air traffic flow. Different mathematical models for optimization strate- 
gies have been elaborated. One approach is a Eulerian model advocated by Menon et al., 
see [MSB04J, where the airspace is divided into line elements corresponding to portions of 
airways on which the density of aircraft can be described as a function of time and of the 
coordinate along the line. This approach focuses on the conservation of aircraft on the line 
elements and uses partial differential equations to describe the time evolution of the process. 
The equations used in this model also appear naturally in highway traffic and were introduced 
by Lighthill-Whitham [LW56] and Richards [Ric56j. This Eulerian network model of air traf- 
fic flow has been considered in several works, e.g. |MSB04| . |BRT04al |BRT04b, BRT06J, 
|SM05], [RSWB06] . |SSBD7| . We also refer to the monograph by M. Garavello and B. Pic- 
coli jGP06j where networks of interconnected roads are modeled and studied, and where the 
Lighthill-Whitham- Richards model is considered on network structures (junctions). 

We use a simplified linear Eulerian network approach. This fits into our scenario since the 
traffic flow is considered as a transport process of aircraft along the edges of a directed graph 
with boundary conditions in the vertices. In this context, the vertices correspond to different 
destinations (or airports) or to bifurcation points of routes in the sky, and the edges model 
the given connections between them (the above mentioned line elements). 

5.2. The allocation matrix. In the literature (e.g., [BRT06J, [SSB07J), the transport pro- 
cesses are usually studied only on an isolated junction of the network. An example is given 
in Figure [3j showing a junction with two incoming edges ei,e2 (called links in [SSB07J) and 
three outgoing edges e3, e^, e$. 



The relation between the incoming and outgoing air traffic flow at a junction is prescribed 
by the so-called junction allocation matrix M(t) = (m,j(t)) for 1 < i < p andp+1 < j < p+q, 
where < m,ij(t) < 1 denotes the proportion of aircraft from incoming link i going to the 
outgoing link j at time t, and 




v, 



Figure 3. A junction with two incoming and three outgoing edges. 



(6) 




j=p+i 
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is required for all 1 < i < p and t > (see e.g. |SSB07| Section 2.3] for this definition). 

As an example, for the allocation matrix M(to) of the junction in Figure [3] at time to, we 
choose 



M(to) 



e 3 e 4 e 5 

\- 4* "i - 

'1/2 1/3 l/6\ <- ei 
1/4 3/4y <- e 2 

This means that at time to, half of the airplanes arriving from edge ei continue their way on 
edge e3, one third of them chooses e 4 and the remaining ones travel to edge es, whereas none 
of the airplanes coming from edge e 2 go in the direction of e3, but one forth of them to e 4 
and the remaining three- fourths to es. Note that Q holds. 

5.3. Our setting and results. We now consider a strongly connected directed graph G 
consisting of n vertices and m edges, and some air traffic flow on it, which, according to the 
linear Eulerian model, can be considered as a transport process. The boundary conditions of 
this process are given in the (transposed) network allocation matrix which we now define on 
the whole network (not only on a single junction) as 

M(t) := (m kl (t)) for k, I € {1, . . . , m}, t € R, 

where < rrikiit) < 1 denotes the proportion of aircrafts arriving from edge e; leaving into 
edge efc at time t. We imply that the flow only takes place on the edges of the network G 
which is given by the (transposed) adjacency matrix of the line graph B = (bki) m xm and we 
set 

m k i(t) = if b ki = 0. 

We further require that 

771 

(7) y^m^ft) = 1 for all I € {1, . . . ,m} and t > 0. 

k=l 



This assumption corresponds to 2.1[ 1 ) (the index / now runs on the edges of the graph) and 



makes the allocation matrix M(i) column stochastic. 

Since every edge of the network only has one end point and one st artin g point, the transpose 



M T (t) of every junction allocation matrix M(t) given in Section 5.2 is a submatrix of the 
bigger matrix M(t). Hence our network allocation matrix M.(t) contains all the information 
stored in the separate junction allocation matrices and ([7| corresponds to the equations ^ 
in every junction. 



We now model the transport process in the network as in Section 2.2 and obtain the 
following air traffic flow problem. 

f t Uj(x,t) = J^Uj(x,t),xe(0,l),t>s, 

(ATF) < Uj (s,0) = (0,1), 

Uj (1, t) = Y^k=i ™>jk(t)uk (0, t),t>0, 

for j = l,...,m. Observe that our non-autonomous boundary conditions together with Q 
imply the Kirchhoff law 2.2 (3j) in the vertices. 



We now proceed as in Section |4j Problems (F) and (ATF) have different solutions since 
their boundary conditions differ. The conditions in (ATF) contain more information and 
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are more demanding. We could look at the problem (ATF) as a subproblem of (F) which 
can also be implemented by constructing a larger graph (the precise implementation is done 
in |Dor08al Section 7.3.2]). In this way we would easily obtain well-posedness, however the 
formulae for the solutions and period would relate to the artificially created larger network 
and it would be difficult to relate it to the original problem. Therefore we rather repeat 
the steps taken for (F) in Section |4j now for the problem (ATF) instead. We will see that 
the main difference is that the weighted adjacency matrix M w (t) is replaced by the network 
allocation matrix M(t). These matrices are different but share many important properties 
(such as positivity, irreducibility, etc.). 

We will assume that the entries of M(i) vary in an absolutely continuous and periodic way. 
This assumption is natural if we think of periodically changing flight schedules (day-night 
rhythms, daily or weekly periods), and without loss of generality we may assume that the 
period is 1. We also assume that the network Gt remains strongly connected at all times i, 
even if some edges of G might not carry any flow at certain times. 

The state space of this system can be modeled by the Banach space X = L 1 ([0, 1], C m ) , 
and the transport process can then be described via the difference operators 

A M (t) :D(A M (t))^X 

defined by 

A M (t)f := f with domain D (A M (t)) := {/ e W 1 ' 1 ([0, l),C m ) \ f(l) = M(t)/(0)} 
for t € R. As in Proposition [5j we can see that the non-autonomous abstract Cauchy problem 

u (t) = A M (t)u(t), t>s, 
u(s) =f s eX. 

corresponds to the transport problem (ATF). Applying our results from Section [3j we obtain 
the following well-posedness result together with a description of the asymptotic shape of the 
solutions. 

Theorem 9. Under the above assumptions, the non- autonomous transport problem (ATF) 
is well-posed. Its unique classical solution is given by the flow evolution family as 

u(t) = U(t, s)f s = M k (t + x)f s (x + t-s-k), 

where f s is the initial distribution of aircraft flow, t > s, x £ [0,l],k<x-\-t — s<k + l, and 
k£N . 

The flow evolution family converges uniformly to a periodic positive group in the sense of 
Proposition [^] with period 

t = lcm {gcd{/ | ej 1 , . . . ej t form a cycle in Gt} \ t E K} . 
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